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Abstract — In CDMA systems, the received user powers vary 
due to moving distance of users. Tlius, the CDMA receivers 
consist of two stages. The first stage is the power estimator and 
the second one is a Multi-User Detector (MUD). Conventional 
methods for estimating the user powers are suitable for under- 
or fully-loaded cases (when the number of users is less than or 
equal to the spreading gain). These methods fail to work for over- 
loaded CDMA systems because of high interference among the 
users. Since the bandwidth is becoming more and more valuable, 
it is worth considering overloaded CDMA systems. In this 
paper, an optimum user power estimation for over-loaded CDMA 
systems with Gaussian inputs is proposed. We also introduce a 
suboptimum method with lower complexity whose performance 
is very close to the optimum one. We shall show that the proposed 
methods work for highly over-loaded systems (up to m (m + 1) /2 
users for a system with only m chips). The performance of the 
proposed methods is demonstrated by simulations. In addition, 
a class of signature sets is proposed that seems to be optimum 
from a power estimation point of view. Additionally, an iterative 
estimation for binary input CDMA systems is proposed which 
works more accurately than the optimal Gaussian input method. 



I. Introduction 

IN a CDMA system, every user is assigned a signature 
vector that is used for transmitting data through a common 
channel. In the channel, all the transmitted vectors are added 
up and the noisy sum reaches the receiver end. 

Because of the variable distance of users from the base sta- 
tion, the received powers of the users differ from their desired 
values. This phenomenon is called the near-far effect; power 
control techniques are used to combat this effect. However, 
because these techniques are imperfect, the receiver needs to 
estimate the received power for each user for proper decoding. 
This problem becomes more critical when the system works 
in an overloaded environment in which the number of users 
is greater than the spreading gain. Thus, the general structure 
for a CDMA receiver consists of two parts, the first one is 
a power estimation module and the second one is an MUD 
module. 

So far, for combating the near-far effect in CDMA, efforts 
have been devoted to finding near-far resistance detectors. 
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MUD rather than single user detection can significantly im- 
prove the performance [1]. However, because optimum de- 
tector is very complex actually NP complete, suboptimum 
detectors are considered. Some of these suboptimum detectors 
are called linear and constrained optimization methods [2]- 
[3]. Subspace-based signature waveform estimation by using 
Wiener [4], [5] and Kalman [6] are other approaches. An- 
other method is based on parametric signal modelling and 
signal spectrum estimation [7] -[8]. Some other strategies for 
combating the near-far effects are the covariance method 
in combination with successive interference cancellation [9], 
blind adaptive interference suppression [10] and isolation 
bit insertion [11]. Additionally, upper and lower bounds for 
near-far resistance of a covariance detector is found in [12]. 
Another approach to mitigate the near-far effect is the channel 
estimation at the receiver Adaptive filter techniques have been 
used to estimate the channel in [13]. However, in overloaded 
systems, the performances of these methods degrade due to 
high multiuser interference. 

Recently, some studies have been performed for CDMA 
systems that do not need a power estimation unit at the 
receiver The discussion in [14] shows that the performance of 
such systems depends not only on the receiver structure but 
also on the signature sequences that are allocated to each user. 
Signature matrices that guarantee errorless communication in 
overloaded systems, when the channel has near-far effects 
in the absence of noise, are studied in the same paper In 
[15], assuming that the near-far effects of the users have a 
Gaussian distribution, the authors have determined some upper 
and lower bounds for the sum channel capacity of the large 
scale binary CDMA systems. 

Nonetheless, the model that has been used for the channels 
in [14] and [15] are much worse than what occurs in practice. 
They assume that the stochastic process of the power change 
of each user is a white process, and the received value of the 
power of each user is independent of its value in any other 
time index. Yet, in practical situations, the fact is that the user 
powers change slowly (in comparison with the data rate); thus 
we can assume that it has a piecewise-constant behavior 

This paper proposes the Maximum Liklihood (ML) esti- 
mation for the user powers even for very highly overloaded 
CDMA systems when the user data and the channel noise have 
Gaussian distribution which achives the channel capacity. The 
method achieves the exact values of the powers asymptotically. 
The only information that is used by the receiver is the 
signatures of the users and the covariance matrix of the channel 
noise. It will be shown that given a CDMA system with chip 
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rate m and no more than m{in + l) /2 users, there exists 
a signature matrix for which the ML estimation has unique 
solution and yields perfect power estimation asymptotically. 
For the systems that have more than m (m + 1) /2 users, the 
ML estimation has no longer a unique solution. 

In addition, two suboptimum estimators are proposed. The 
first one is a suboptimum version of the proposed ML esti- 
mation that has lower computational cost with similar perfor- 
mance. Additionally, this study will introduce a characteristic 
about how a signature matrix is good from power estimation 
perspective and a new class of codes will be proposed that 
seems to be optimum with respect to this characteristic. 

The second method is an iterative power estimation in case 
of binary input vectors. This iterative method works much 
more accurately with just a few number of received vectors. It 
will be proved that for a noiseless channel with user powers 
constrained to be in a predefined interval, this method can find 
the powers exactly. Simulation results show that for high noise 
variance and wide range of powers, the estimator remains 
accurate. 

The performances of the proposed methods are simulated 
for different situations when the channel noise variance is 
known or unknown at the receiver The ability of tracking the 
changes of the powers is also simulated. Moreover, a CDMA 
system is simulated that uses binary signatures in conjunction 
with BPSK modulation utilizing the proposed methods for 
power estimation and is also compared to the system with 
perfect or no power control. The simualtion results show 
excellent performances for these methods. 

In the next section, the channel model is presented. In 
Section Uni the ML power estimation is introduced. In Section 
IIVI a suboptimum estimation is proposed. We also introduced 
a class of accurate signature matrices for estimating user 
powers. The iterative estimation is discussed in section |V] 
These two methods are simulated in a practical system using a 
BPSK modulation in section lVTl Conclusions and future works 
are in section IVIII 

II. Channel Model 

In a synchronous CDMA system, each user is assigned 
an TO-chip signature vector for transmitting its data through 
the channel. Each user multiplies its data by its signature 
vector before transmission. The channel noise is added to 
the transmitted vectors with different powers at the receiver. 
Assuming a synchronous CDMA system with n users such 
that the received power of the z*'' user is pi, i = 1, ... ,n, the 
channel can be modeled as 



Y = SP^/^X + N 



(1) 



where S is the mxn signature matrix in which the i*'^ column 
is the signature of user i, P = diag (pi, . . . ,p„), X is the 
user data vector, N is the channel noise vector, and Y is the 
received vector Notice that it is assumed that the columns of 
the signature matrix S are normalized and the entries of X 
have unity power 

Because of the variable distance of the users from the base 
station, the values of the user powers p^'s change and thus the 



matrix P is unknown at the receiver Therefore, at the receiver 
we need to estimate the user powers for proper extracting data 
X from the received vector Y. 

In practical situations, the data rate is high enough such that 
we can reasonably assume that P is constant over a number 
of received vectors, say L. In other words, the variation of the 
matrix P versus time is piecewise constant. 

Power estimation is not difficult in under- or fully-loaded 
systems (n < m) since orthogonal signatures that provide an 
independent channel for each user can be used. However, in 
over-loaded systems (n > to), where such codes do not exist, 
the problem becomes much more complicated. Nevertheless, 
over-loaded systems are preferred because of their bandwidth 
efficiencies. The following section proposes the ML estimation 
of the matrix P for a system with any given number of users 
with Gaussian distribution which achives the channel capacity. 

III. ML Power Estimation 

In this section using ([T]), we wish to find the maxi- 
mum likelihood estimation of the user powers matrix, P = 
diag {pi, . . . ,pn), from the observed vectors Yi, . . . ,Yl. We 
assume that in these L observations the user powers are 
constant. Thus, we have 



ML 



arg max / ( Yi , 



, YA F) 



(2) 



where Pml is the ML-estimation of the user powers matrix 
and / is the Probability Density Function (PDF) of the 
observed vectors. Hereafter in this section, we assume that the 
data vector has Gaussian distribution with covariance matrix 
E{XX^} = Im and the noise vector N has distribution 
AA(^,S). 

Since F/s for i ~ 1,...,L are i.i.d. Gaussian random 
vectors with mean fi and covariance matrix of Mp = SPS"^ + 
S, we have 

argmax []/ P) = argmax In ([] / (^'1 P) ) 

i=l \i=l / 

In /mi P) 

(3) 



ML 



L 

= arg max ' 



: arg m^x 



lnJ(2^)'"|det (Mp) 



Thus, we deduce that. 



ML 



arg mm ( | - f^f Mp^ (F, - 

^ i=i 

+ In |det (Mp) 



(4) 



Notice that the estimation of the user powers is possible if 
the solution of the above equation is unique. The following 
theorem expresses conditions under which the above 
minimization has a unique solution. 
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Theorem 1: Let S be an xn matrix whose rows are 

the entry-by-entry multipHcation of the rows of the signature 
matrix S. Then, the minimization (|4|i has a unique solution if 
and only if 



rank (^S j = n (5) 
We prove the necessary condition of Theorem [T] and the 
converse will be proved for sufficiently large i's. 

Proof of necessary condition: Suppose that n > rank 
Then, there exists an n x 1 non-zero vector B such that 
SB = 0. It yields that SBS^ = where B IS an n X n 
diagonal matrix that has the entries of B_ as its diagonal entries. 
Now, since for any diagonal matrix P, Mp = SPS^ + S = 
S (P + B) + S = Mp+B, the minimization (|4| has no 
unique solution. 

Proof of the converse theorem: We desire to show that for 
sufficiently large L, the second derivative of (|4]i is positive 
definite and thus it is a convex function and has a unique 
minimum. Let 

1 ^ 

* (P) = 7 ^ {Y^ - tif Mpi (r, - /i) + In |det (Mp)| 



Then we have. 



(6) 



^ (A, P) = Dp* (A) = tr (Mp'^SAS^) 



1 

- ^ (r, - ^if Mpi (SAS^) Mpi [Y, - /i) 



= tr(MpiSAS^) 



= tr(S^MpiSA) 



tr I^S^Mpi Yl - - '^)'^) MpiSA^ 

= tr (S^Mpi (I - WMpi) SA) (7) 

where _Dp$ (■) denotes the derivative of $ at point P, W 
is the estimation of the covariance matrix Mp from the re- 
ceived vectors, i.e. W — j- X^i^Li (Yi — m) {Yi — Getting 
another derivative, we obtain 

L>P*(A,B) = 

tr(S^MpiSAS^Mpi(2W-Mp)MpiSB) (8) 

where A and B are diagonal matrices. For the proof of 
convexity, we wish to show that Dp^{A,A) > and 
equality holds only for A = 0. Note that 

Dp^- (A, A) = tr (Mp^AMp^ (2W - Mp) Mp^ a) (9) 

where A = SAS^. We have 

tr (^Mp^AMpi (2W - Mp) Mp^A^ 
= tr (Mp^AMpiAMpi (2W - Mp)) (10) 



Since 

Mp^AMpiAMpi = Mp^/^ (Mp^/^AMp 



-1/2 



'Mp'/' (11) 



— 1/2 

and Mp ' AMp*'" is a symmetric matnx, 
Mp^AMp^AMp' is a positive-semidefinite matrix for 
non-zero A. As L tends to infinity, from law of large 
numbers, 2W — Mp tends to Mp which is a positive-definite 
matrix. Therefore, for large L, according to Lemma [T] 
described below, tr (Mp^AMp^AMp' (2W - Mp)) > 



-1/2 



and equality holds when A = 0. Since rank ySj — n, A = 
is equivalent to A = and thus the proof is complete. ■ 

Lemma 1: If V and W are positive-semidefinite matrices 
such that ^ 0, then tr (VW) > 0. 

Proof: tr(VW) = tr (vVSwVi/^). Since W is a 
positive-semidefinite matrix, V^/^WV^^^ is also a positive- 
semidefinite matrix. Thus, tr(VW) = tr (Vi/^wVi/^) > q 
and equahty holds only when V^/^WV^/^ = 0.1 

Theorem [T] shows that for a given signature matrix S, the 
rank {s^ is an upper bound for the number of users for which 
the user powers can be estimated from the received vectors. 
Obviously, rank ^s) < I2(!!i±l) the next sections, it will 

be shown that for any to, there exist to x "'(™+^) signature 
matrices S such that its corresponding S is full rank. 

By having a system with chip rate to, it is possible to 
estimate the user powers for up to to (to + 1) /2 users. It is 
interesting that if the number of users is beyond this threshold, 
the user powers cannot be estimated from the received vectors 

Iy's. Also, note that if n < m (to + 1) /2, this estimation is 
asymptotically exact, i.e., if the number of received vectors 
tends to infinity, the estimated powers tend to the exact values 
with probabiliy 1. This fact is very significant because it 
implies that we can estimate user powers for any practical 
overloaded CDMA systems since the number to (to + 1) /2 is 
way beyond number of users for any practical situations. 

In addition, using (|7]) in the proof of Theorem [T] we can 
implement the ML power estimator If P = V ml minimizes 
(IHi, for any diagonal matrix A, we have 

tr (S^M-i (I - WM-i) SA) 



0. 



(12) 



It is equivalent to that all diagonal entries of 
S^Mp' (I - WMpi) S are zero, i.e., 

[S^Mpi(l-WMpi)S]^^ = 0, (13) 

for i = We can solve this system of non-linear 

equations to obtain the ML estimation of the user powers. 

Fig. [U is the simulation of a system with 12 users and 8 
chips with binary WBE signatures [17]. We have used the 
Newton method for solving the above system of equations. 
The simulation is for AWGN channels with two specific values 
of Eb/No. 

In the next section, we propose a suboptimum estimator that 
has a lower computational complexity. Also, we introduce a 
family of signature matrices which are good from the power 
estimation point of view. 
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Fig. 1. The ML estimated power vs. the number of received samples for 
Et/No = 6dB and 16dB. 



Fig. 2. The ML and suboptimum estimation of the powers vs. the number of 
received vectors L for a system with n = 12 and m = 8 that uses a binary 
WBE signature matrix. 



IV. SuBOPTiMAL Power Estimation 

In this section, a low complexity method for estimating pow- 
ers is introduced. Simulation results show that the performance 
of the proposed low complexity estimator is very close to the 
ML estimation introduced in the previous section. 

A. The Proposed Method 

Since W and Mp are respectively the empirical and the 
exact covariance matrices of Y, it is expected that W and 
Mp are approximately the same. Notice that, using the law 
of large numbers, W tends to Mp as L tends to infinity. 
Hence, we try to find the diagonal matrix P, such that 

W = SPS^ + s 



for unknowns, pi , 



We rewrite (fl4] i in the form 



K ^SP 



(14) 



(15) 



where K is an 



X 1 vector which contains the entries 



m(?n+l) 



X n matrix 



of the upper triangle of W — S, S is an 
defined in Theorem [T] and P ^ [pi, ■ ■ ■ ,Pn] In general, we 
consider the least square solution of (fTsl l for P. The same as 
the result of Theorem [T] the least square solution of (ITSi is 
unique if and only if the rank (^S^ — n. 

There is another justification for choosing powers stated in 
the previous paragraph. In the preceding section, we showed 
that the ML estimation of the matrix P should satisfy ( fT2b . 
Here P is the annihilating term for I — WMp^ in ( fT2b . 
Interestingly, simulation results show that the performance of 
this simpler estimator is as good as the ML one. Fig. |2] 

In practical situations, we face two other problems. Firstly, 
the covariance matrix of the channel noise is not always known 
at the receiver end. Secondly, the user powers change by time. 
Using the suboptimum estimator proposed in this section, we 
devise solutions for these problems. 

If S is unknown but the channel noise is white, S is a 
diagonal matrix. In this case, from (fT5t we omit equations 
that correspond to the diagonal entries of S. Thus, we arrive 
at a system of '"("'~^) equations. 



where S is an 



X n matrix wherein each of the rows 



is the entrywise multiplication of different rows of S (the 
entrywise multiplication of a row by itself does not appear 
in S) and is the "'(™~^) x 1 vector that contains all 
entries of W that are above its main diagonal. Again, in 
general, we choose the least square solution of (fT&i as the 
estimated powers. This solution is unique if and only if the 
rank ^S^ = n. It means that for a system with white noises 
and unknown variances, the maximum number of users cannot 
be more than m (m — 1) /2. 

To add the ability of tracking the user power variations, we 
use another window for calculating the empirical covariance 
matrix of Y. By assigning 



W 



1 



L 



(17) 



K ^SP 



(16) 



where < a < 1, we decrease the effect of the old samples 
of the received vectors. Figs. |3] and |4] show that by this 
modification, the estimator tracks the power changes rapidly. 
While it almost perfectly tracks the sinusoidal fluctuation, it 
is interesting to see that it can also track abrupt changes. 

In the next subsection, we present the signature matrix S 
whose corresponding S and S have rank n, respectively. More- 
over, the proposed matrices may lead to maximum robustness 
against errors in W. 



B. Signature Matrix Design for Estimating User Powers 

It is known that sensitivity of P to the changes in Kjs equal 
to the reciprocal of the minimum singular value of S. Thus, 
for the proposed estimation method to be more robust against 
the error in estimating K_, a matrix S with larger minimum 
singular value is needed. More important, if the rank of matrix 
S is less than n, its minimum singular value is and thus ( fTSl ) 
lacks a unique least square solution. The same statements hold 
by substituting S for S and for K_ in the system where the 
noise variance is not known. 
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Fig. 3. The estimated power vs. the number of received samples for the 
proposed suboptimum estimator for n = 12, m = 8 and a binary WBE 
signature matrix; the user powers are changing with a sinusoidal shape. 
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Fig. 4. The estimated power vs. the number of received samples for the 
proposed suboptimum estimator for n = 12, m = 8 and a binary WBE 
signature matrix; the user powers is doubled abmptly. 



In the following theorem, it will be shown that for any 
m and n < m_{m+ 1) /2, there exists S^xn such that its 
corresponding S m(m+i) has rank n. 

Theorem 2: Suppose S^xn is a matrix that its coresspond- 
ing S ^(m+i) has rank n. Also, for any fc < m + 1, let 



S(rn+l)x(n+i;) 

Then S(^+i)(m+2) , 



llxn 







ifexfe 

(m+l — k) X k 



(18) 



has rank n + k. 



we find that k 



■ X (n+k) 

Proof : Constructing S (^+ik^+2) ^ 
linearly independent columns are added to the columns of 
Sm(m+i) . Using the assumption that S n.(m+i) has rank n, 

2 ^ ^ 2 ^ 

S (m+i)(m+2) , , , has rank n+k and the proof is complete. ■ 

The following theorem shows that for any m and 
n < TO (m — 1) /2, there exists S„ixn such that its 
corresponding S m(m-i) has rank n. 

2 ^" 

Theorem 3: Suppose Smxn is a matrix that its coresspond- 



ing S ,ri(ni-i) has rank n. Also, for any fc < m, let 



'(m+l)x(n+/£) 



Olxri 
Q 

X r?. 



■-fcxfc 







(rn — k) X fe. 



(19) 



Then S . . , ^ has rank n + k. 

The proof is simillar to that of Theorem |2] 

Corollary 1: Let Sixi = [1]; using Theorem |2l we can 
construct appropriate matrices for power estimation for any 
TO and n < m (to + 1) /2. Also, using Theorem |3] let 
§2x1 [l l] , we can construct suitable matrices for power 
estimation for any to and n < to (m — 1) /2. 

Since multiplying each column of S with a non-zero number 
leaves the rank of S and S unchanged, the columns of the 
matrices introduced in Corollary [Tjcan be normalized to obtain 
acceptable signature matrices. 

An interesting fact about the normalized version of the 
matrices constructed in Theorem |2] is that their S have 
large minimum singular value and thus are almost optimum 
for estimating the user powers. For the second set of the 
constructed matrices, in the following theorem we prove that 
all singular values of the S are 0.5. 

Theorem 4: Let Smxn be a signature matrix defined in the 
Theorem [3] All singular values of the S are equal to 0.5. 

Proof: From the definition of S, it is deduced that S is 
formed by concatenating a permutation matrix multiplied 
with 0.5 with an all zero matrix. Hence, all the singular 
values are 0.5.H 

It is worth mentioning that by checking more than 50000 
normalized random matrices, we did not find any matrices with 
minimum singular value greater than the ones constructed by 
Theorems |2] and [3] However, the optimality of these matrices 
has not been proven yet. 

The performance of the matrices proposed in Theorems 
and [3] are simulated. The simulations were performed for the 
proposed 8 x 36 and 8 x 28 matrices in situations where the 
noise variance is known and unknown at the receiver end, 
respectively. Although the overloading factor is very severe, 
the powers are estimated accurately by a few number of re- 
ceived vectors. Notice that the conventional power estimation 
methods fail to work even in much lower overloaded systems. 

V. Iterative Power Estimation for Binary Input 
CDMA System 

In pevious sections, we introduced optimum and subopti- 
mum power estimators in the case of Gaussian input vector 
However, in practical situations, binary input systems are 
more favorable. While estimators proposed in Sections HIT] 
and HV] work for binary input, they are not really optimum. 
In this section, iterative power estimation for binary input 
synchronous CDMA systems is developed which is superior 
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Fig. 5. The comparison of the estimated power vs. the number of received 
samples for an 8 X 36 matrix (Theorem |2) and a random matrix of the same 
size. Here we have an AWGN channel with a known noise variance and 
Et/No = 6dB. 
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Fig. 6. The comparison of the estimated power vs. the number of received 
samples for an 8 X 26 matrix (Theorem |3) and a random matrix of the same 
size. Here we have an AWGN channel with an unknown noise variance and 
Ei/No = 6dB. 



to the estimators described previously. Estimations become 
remarkably better when Eij/Nq is increased. 

Indeed it will be proved that with no noise in the channel 
this method can determine the actual user powers with proba- 
bility one provided that the variance of the power fluctuations 
do not exceed a given threshold. However, the simulation 
results show that this method has an acceptable performance 
in noisy channels and when power fluctuations are very large. 

We use an iterative scheme for estimating the powers. 
In each iteration, we use the output power of the previous 
iteration and decode the user data. Then, using these extracted 
values of the user data, we again estimate the user powers. Fig. 
|7] shows this procedure. We denote the decoded data of the i*'* 
received vector and the power matrix in the iteration by 
Xl and P^, respectively. 

According to the supposition that the powers are constant 
in the interval of sending L vectors, we find the least square 
solution of the following linear system of Lm equations and 
n unknowns (p^'s for i — 1, . . . ,n) as the estimation of the 
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Fig. 7. The j* iteration of the iterative power estimation 



1,2, 




(20) 



By Multiplying each column of the signature matrix S by its 
corresponding entry of Xl, we get S^. Thus (l20l i become 



(21) 







"Si- 




Lm X 1 


Si 



Lm X n 



We should of course guarantee two main postulates in the 
above procedure. Firstly, we should make sure that the set of 
linear equations in (l2Tl i has a unique least square solution or 
in other words the matrix has rank n. Secondly, it should 
be proved that the above procedure will converge to the actual 
value of user powers. In the following, we express conditions 
in which the convergence of iterative method is guaranteed. 

Theorem 5: Suppose S contains a row that has no zero and 
Xi's are uniformly chosen from {±1}". Then the probability 
that Simxn — [Sf ■ • ■ S|^]"^ has rank n tends to one as L 
approaches infinity, where S,; is formed by multiplying each 
column of the signature matrix S„ixn by the corresponding 
element of Xi. 

Proof: Suppose that the q*'* row of S has no zero. Form 
an L X n matrix C whose fc*'* row is the (km + q)*^ row of 
S for fc = 0, . . . , (L - 1). Since the rank(S) > rank(C), it is 
sufficient to show that the rank(C) — n with probability one 
for sufficiently large values of L. This matrix can be written 
as: 



C = 



-xf 




Sg,l 






^q.n 



(22) 



where Sg_,; is the i 



th 



entry of the q row of S. Note that the 
second matrix is diagonal with rank n. Thus it suffices to show 
that the rank of the first matrix approches n as L increases. Let 
Pl be the probabilty of having n linearly independent vectors 
among L Ixn random vectors. By partitioning C into [L/n\ 
blocks we have the following inequality: 

PL>l-il-Pn)^''^'"^ (23) 

where P„ is the probability of a random n x n matrix with 
{±1} entries being invertible. The probability P„ is lower 
bounded by [19] 



(24) 



Pn>l-{j+ 0(1) 
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Fig. 8. Estimated user powers using the iterative method with 4 iterations 
for a noiseless channel, rj = 0.5 . 



Thus we have: 



Pl > 1- ( ^ + 



n [L /nj 



(25) 



According to (|25l l tends to unity as L approaches infinity 
and the proof is complete. ■ 

According to [14], for any uniqely decodable signa- 
ture matrix there is an rjsup such that if pi's belong to 
[1 ^ Vsup, 1 + Vsup], then the user data can be extracted with- 
out any error. The following theorem shows that in this situa- 
tion the user powers can be determind exactly with propability 
one. 

Theorem 6: For any uniquely decodable signature matrix 
and the corresponding rjsup explained in [14], if p/s belong 
to [1 — rjsupA + Vsup], then in a noiseless channel we can 
exactly find user powers with probability one as L approaches 
infinity. 

Proof: According to [14], in the first step of the above 
procedure, X/'s are exactly the user data which are random. 
Therefore the set of equations in (l2ll are consistent, i.e. it 
has an exact solution. Moreover, Theorem |5] states that as 
L increases, the probability of having a unique least square 
solution in (ISTT i approaches unity. Hence the resultant powers 
after the first step are the desired results. ■ 

In the following, we have simulated the above method. In 
the simulations, the signature matrix is an 8 x 13 COW matrix 
proposed in [18] and the initial value for the powers in the 
first iteration is unity. Fig. |8] shows the result for a noiseless 
channel. Although according to [14], the corresponding ris^p 
is equal to 0.23, we simulated the system for 77 = 0.5. Notice 
that although Theorem |6] guarantees perfect estimation for 77 < 
0.23, Fig. [8] depicts that even for higher fluctuations in the user 
powers, the performance of the method is impressive and the 
powers are derived after 4 iterations. 

As previously stated, the constraints introduced in Theorem 
|6]is not the only situation that this method works. Simulation 
results prove that this method is not only useful for noisy 
channels but it is also applicable for high fluctuations of the 



Fig. 9. Estimated user power using the iterative method with 4 iterations. 
In this case user power is a Gaussian random variable with mean 25 and 
standai'd deviation 10. This simulation is done for two different values of 
Et/No, 6dB and 16dB. 
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Fig. 10. Estimated user power using the iterative method with 4 iterations. 
In this case user powers are doubled abruptly and Ed/No = 6dB. 



user powers. It is expected that the effect of noise can be 
somehow compensated by finding the least square solution 
as shown in Fig. |9] In this case, the user power is a Gaussian 
random variable with mean 25 and standard deviation 10. This 
simulation is performed under two different values of Eh/No, 
6dB and 16dB. Interestingly, the powers are obtaind using only 
30 received vectors. 

Simulation results show that the variance of the user powers 
is critical and their means have nothing to do with the 
estimaion procedure. 

Figs. [TOl and [TTI show the ability of this method in tracking 
the user power changes. These curves are obtained with only 4 
iterations and E^/No — 6dB. We have used a window which 
keeps the last 40 received vectors for estimation. As you can 
see this method estimate powers more rapidly and accurately 
in comparison with suboptimum method. 

VI. Bit-Error-Rate Performance 

In this section binary CDMA systems that use the subopti- 
mum and iterative power estimations have been simulated. A 
64 X 104 COW matrix [18] as the signature matrix is used. 



8 



3.5 




Estimeted Powei 



Actual Power 

qI , , , , 

200 400 600 800 1000 

Sample Index 

Fig. 11. E.stimated user power using the iterative method with 4 iterations. 
In this case user powers have sinusoidal shape and Ei,/No = 6dB. 

The simulation is performed for four different cases. In the 
first case, user powers are known at the receiver end, i.e., 
perfect power control. In other cases the performance of the 
system is simulated for suboptimum and iterative estimators 
and they are compared to the no power control case. 

These simualations are performed for two different cases. 
In the first one, user powers are modeled as a sinusoidal 
changing variable in the range of [0.5, 1.5] with period of 
20000 samples. In the second case, user powers have piecewise 
constant variation with length 10000 samples and abrupt 
changes to randomly chosen amplitudes from the [0.5, 1.5]. In 
fact, these fluctuations for user powers are much faster than 
what appears in practice. 

In [18], a very fast ML decoding scheme has been proposed 
for COW codes. This method can be implemented with a 
similar number of computations for the case of the near-far 
effect. In order to make the decoder of [18] suitable for the 
proposed model, it is necessary to make a minor adjustment. 
After estimating the powers pi, . . . , pn the decoding algorithm 
has the following two steps: 

. Step 1: Let S„ixn = [A™xm|B„x(ri-m)] where A is an 

invertible matrix. 
. Step 2: Assume X = [Xf Xj]^ e {±1}" where Xi 

and X2 are to x 1 and (n — m) x 1, respectively. Find 

X2 that minimizes 

\\A-^Y - A-IBP2X2 

- Pisign(A-ir - A-iBP2^2) II (26) 

and set 

Xi = sign(A-iy - A-^BP2^2) (27) 

where 

Pi =diag(pi,...,p„) (28) 

and 

P2 = diag {pn~m+l, ■ ■ ■ ,Pn) ■ (29) 

Figs. [12] and [T3] show the simulation results for sinusoidal 
and stepwise fluctuations, respectively. We use only 4 itera- 
tions in these figures for the iterative method. 
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Fig. 12. Comparison of BER vs. Et/No for the BPSK modulation CDMA 
systems using the suboptimum and iterative power estimations when user 
powers have sinusoidal fluctuations; the extreme cases of perfect and no power 
controls are also depicted. This CDMA system consists of 64 chips and 104 
users. 
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Fig. 13. Comparison of BER vs. Eb/NO for the BPSK modulation CDMA 
systems using the suboptimum and iterative power estimations when user 
powers have stepwise fluctuations; the extreme cases of perfect and no power 
controls are also depicted. This CDMA system consists of 64 chips and 104 
users. 

Note that both methods saturate for high values of Eb/NQ. 
This is because when the power changes, it takes some time 
for both methods to track the power variations. In this period, 
we always have some errors whether there is noise or not. 

VII. Conclusion and Future Studies 

In this paper, the optimum power estimation (i.e. ML 
estimator) for a given model of stochastic process of user 
powers has been discussed. We have also introduced a lower 
complexity power estimation method which can perfectly 
estimate the user powers in a highly overloaded CDMA 
system. In fact simulation results have shown that the ML 
estimation and its simplified low complex version have similar 
performances. This paper also introduced a new characteristic 
for the signature matrices that shows their suitability for 
estimating the user powers. It was shown that for a system 
with TO chips, the user powers are not attainable by the ML 
method unless the number of users is less than m,{m + l)/2. 
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This characteristic is independent of the previously known 
characteristics such as WBE [16] or COW [18] and thus it 
should be incorporated seprately when designing the signature 
sequences for a CDMA system. While these methods do not 
use the knowledge of the input alphabet of the users, we 
have proposed an iterative method that works for binary input 
vectors and its performance is much better than the previous 
method. 

Finding the optimum signature sequences according to the 
proposed characteristic for every chip rate m and the number 
of users n is an interesting topic for future research. By 
exploiting this knowledge, it would be possible to achive more 
efficient power estimation methods for binary input systems. 
The extension of this work to asynchronous CDMA systems 
would be another worthwhile research. 
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